Single-file diffusion with non-thermal initial conditions 
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Single-file diffusion is a theoretically challenging many-body problem where the calculation of 
even the simplest observables, e.g. mean square displacement, for a tracer particle requires a heavy 
mathematical machinery. There is therefore a need for simple approaches which predict qualitatively 
correct behaviours. Here we put forward one such method which we use to investigate the influence 
of non-thermal initial conditions on the dynamics of a tracer particle. With our new approach 
we reproduce, up to scaling, several known asymptotic results for the tracer particle mean square 
displacement. 



I. INTRODUCTION 

The dynamics of diffusing particles in one dimension 
which are unable to pass each other, often referred to as 
single-file diffusion, has attracted great theoretical inter- 
est for at least 50 years, see [THI1|. More recent progress 
in the field addresses single-file diffusion in a potential 
[15j and that of particles with different diffusion con- 
stants [16]. One of the main theoretical predictions is 
that a tagged (or a tracer) particle, explores the sys- 
tem subdiffusively even though the collective behaviour 
is identical to that of non-interacting particles [121 HZ] ■ 
The subdiffusive behaviour can be understood through 
the tracer particle's velocity-velocity correlation function 

3/2 

which is negative and decays as tj^ ' with delay time Id 
[4]; this results translates into the classical square root 
law for the tagged particle's mean square discplacement 
S(€) ex i^/^ where t is time. 

Nowadays single-file diffusion also finds experimental 
realizations. A few examples are diffusion in colloidal sys- 
tems [IH1II9], transport in microporous materials [20ll21j. 
and permeability of potassium ions in nerve fibres |22j . 
Also, exclusion effects have been proposed to be of im- 
portance for transcription factor (DNA binding proteins) 
dynamics [j23 . 

In order to calculate dynamical properties of a tagged 
particle in a single-file system one needs to deal with 
a theoretically challenging many-body problem. This 
typically requires rather elaborate mathematics (e.g. 
[TJ m [m [131 [15]) which may cloud fundamental under- 
standing. There is therefore a need for simplified meth- 
ods that captures correct qualitative behaviour and at 
the same time provide new insights to the problem. A 
few such papers already exist [51 E] for thermal ini- 
tial conditions and identical particles. To the best of 
our knowledge there has not been a similar progress 
on simple approaches for non-thermal initial conditions. 
Some aspects have, however, been addressed [MH^ us- 



ing rather involved techniques. The main goal of this 
paper is to put forward a new but simple method to 
calculate tagged particle properties, in particular S(t), 
under non-thermal conditions. Our approach relies on 
the assumption that we can make a quasi-static approxi- 
mation for the velocity-velocity autocorrelation function 
and set it equal to the correlation function at equilib- 
rium. Based on this we are able to reproduce known 
scaling behaviours pertaining to different types of initial 
conditions. 



II. THE QUASI-EQUILIBRIUM APPROACH 

Imagine a single-file system extending to infinity in 
both directions. The particles are identical, point-like 
and undergo Brownian motion in between hard-core 
repulsive interactions. Now we tag one particle, la- 
bel it "0", and place it at time i = in the ori- 
gin, that is, X[){t = 0) = 0. The remaining particles 
are initially distributed symmetrically around xq = 0, 
...a:_2, a;_i, Xo, cci, a;2, ..., according to some non-thermal 
distribution pq(x). We then ask: what is the mean 
square displacement, S(t), of particle 0? Here we define 
5(i) = {xo{t)'^) where (...) denotes ensemble average. 

Our expression for S(t) is based on the following quasi- 
equilibrium argument. For an equilibrated single-file it 
is well known that the long time limit of the velocity- 
velocity correlation function for a tagged particle is [4] 



{v{t)vit - to)) 
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(1) 



where v(t) — dx{t)/dt and g is the constant average 
particle concentration (average inverse inter-particle dis- 
tance), and D is the single-particle diffusion constant (as- 
sumed equal for all particles). Now we assume that Eq. 
(n]) holds for all, even non-thermal, initial conditions if we 



replace g with the time-dependent concentration p(x, t) 
in the vicinity of the tagged particle, i.e., 



Q —!• p{x — 0, t) 



(2) 



This is a quasi-equilibrium approach which we as- 
sume to hold since concentration relaxations propa- 
gate faster than the tracer particle [23] • Integrating 
Eq. ([I]) twice with respect to time yields S{t) — 

Jldt' Jldt"{v{t')v{t"))=2!ldt' f^ duHt')v{f-u)). If 
we furthermore assume the tagged particle behaves elasti- 
cally over time-scales longer than the velocity relaxation, 
we have /„ du{v{t)v{t — u)) =0 [30] ■ This means that 

we can write S{t) = —2 J^^ dt' J^, du{v{t')v{t' —u)) , which 
together with Eq. (IT]) and g = p{x = 0, t') gives 



(3) 




t'p{x^O,t') 



Now, since the macroscopic particle concentration in a 
single-file system and in a system of non-interacting par- 
ticles are indistinguishable [T7] we can calculate p{x = 
0, t) by simply solving the diffusion equation in one di- 
mension, and subsequently setting a; = in the asso- 
ciated solution. In effect this amounts to convoluting 
^0(2^0) with a Gaussian propagator |31] 



p{x = G,t) 
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(4) 



Equation (pi [together with Eq. (B] constitutes our main 
result. It is straightforward to see that for uniform initial 
densities those equations lead to classical single-file result 
for the mean square displacement 



g y TT 



(5) 



(see e.g. [T]) as it should. 



III. SPECIAL CASES 



In this section we will investigate four different types 
of poixo) and show that our simple model gives correct 
results up to scaling in t. Our results are summarised 
briefly in Tab. |l]and depicted in Fig. [ij 
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TABLE I: Summary of results (r = A^/D). 



leads to 



pix = 0,t)^^e'/-eTfc[^t/T) (7) 

which together with Eq. ([3| gives S{t) predicted by our 
model. Here we introduced 



A2 

Id 



(8) 



Figure [T] (solid blue line) shows the result of a numerical 
integration of Eq. ^ [using Eq. ([7|]. We see clearly that 
S(t) exhibits two regions with different dynamics with 
crossover time r. The limiting behaviours can be found 
analytically if we expand Eq. nh in the short and long 
time limit: for short times {t <^ r) we get p{x = 0, t) ~ 
N/{2A), whereas for long times (t ^ r) p{x = 0,t) ~ 
N/VAttDI. The resulting mean square displacement in 



the two limits reads 







(9) 



This means that for short times the system is crowded 
enough such that we recover the known scaling with time 
in Eq. ^ . For long times the tagged particle follows the 
system's center-of-mass motion yielding S{t) ^ Dt/N. 
For this limit we point out that C. Aslangul [27] obtained 
irDt/N using a more elaborate approach. 



B. Gaussian distribution 



A. Exponential distribution 

Here we consider the case where the initial particle 
density decays exponentially from the origin 
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where N is the number of particles and A is the char- 
acteristic decay length. Using this expression in Eq. ^ 



In this subsection we take the particles to be initially 
distributed as a Gaussian around the origin 
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Putting this expression into Eq. (B gives 

N N I I 

p(a; = 0,t) = - 
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(10) 



(11) 



with T given by Eq. ^. The resulting S{t) from Eq. ^ 
is shown in Fig. [I] (solid red line) where we see, again, 
two dynamically different regions with crossover time r. 
A similar analysis as in the previous subsection shows 
that 



S{t) 



2Dt/N ex i, t > r 



(12) 



where we used that p{x — Q,t) ~ N/{A^/2n) for short 
times and p{x = 0, i) ~ N/^AnDt for long times. The 
crossover time as well as the limiting behaviours of S{t) 
are up to scaling in agreement with (24) . 
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C. Equidistant particles 

Here we address the situation where the particles are 
initially distributed on the line with the same distance A 
apart from each other: 



10" 



Po{xo) = ^ S{xo - nA) 



(13) 



where S{z) denotes the Dirac delta function. This case is 
different in one important aspect compared to the other 
cases in this paper: the concentration p{x, t) in the long 
time limit approaches the (finite) constant 1/A whereas 
in the other cases considered here p{x, t) goes to zero. 
This means that we expect S{t) ~ ^flji for long times 
for the present case. If we carry out the calculation of 
p{x — 0, t) as before we find 



p{x^Q,t) = 



iMVW^ 



A ^47ri/T 
where 'd3{z) is the elliptic Jacobi-Theta function [28 



(14) 
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The behaviour of S{t) is depicted in Fig. IT] (solid green 
line) where two regimes are visible once more. Using that 
p{x = 0, i) ~ l/\/47rZ?t for short times and p{x = 0, t) ~ 
1/A in the long time limit, we obtain the short and long 
time asymptotics: 




t> T 



(16) 



The long time limit is in agreement with [26] up to scal- 
ing. However, the prefactor is a factor v2 too large com- 
pared to the correct result. Our simple model is therefore 
unable to account for numerical factors but gives the cor- 
rect scaling behaviour in all cases we examined here. 

One reason to why our method is unable to produce 
proper prefactors is that the density p{x^ t) does not con- 
tain all necessary information to predict the motion of 



10" 

t/T 



10' 



FIG. 1: Mean square displacement versus time for different 
initial distributions po{xo) (r = A^/D). The curves are ob- 
tained from a numerical integration of Eq. ([3| 



the tagged particle. One way of seeing this is from the 
following simple example. Imagine that all particles are 
placed equidistantly with inter-particle distance A. Now 
we draw a random number between zero and A from 
a uniform distribution and add the same random num- 
ber to the position of all particles. If we average over 
many such subsystems we would obtain a new system 
with uniform density 1/A just as in thermal equilibrium 
(each single-equidistant case has of course Dirac delta- 
peaked density). However, since each system is just the 
equidistant case shifted by a constant, S(t) for the aver- 
aged system will still be a factor ^/2 different compared 
to the case of thermal initial conditions. This illustrates 
that two cases may have the same uniform density but 
yet yield different prefactors to the \/i behaviour. Thus, 
the density p{x,t) is by itself not sufficient to predict 
exact numerical prefactors. Another way of seeing that 
p{x, t) does not hold all necessary information is to con- 
sider S{t) in thermal equilibrium. For this case Alexan- 
der & Pincus [9] argued that S{t) can be related to the 
density-density correlation function (dynamic structure 
factor) rather than the density itself. Their relation is 
another indication that only the density is not sufficient 
to predict the full S{t) but that one also need to consider 
density fluctuations at different locations. 



D. Power-law distribution 

Here we investigate the scale-free initial density 

Po(a;o) = ^ (^) , ? < |xo| < oo, < /3 < 1 (17) 
where we take the cut-off ^ — >■ 0. The power-law form 



(17) yields standard single-file dynamics [Eq. (|5|] for 



/3 = 0. If we put Eq. ([T7| into Eq. Q we find 



kH'-^)i^r <-' 



where T{z) is the Gamma function. Using this formula 
in Eq. (|3| leads to a mean square displacement charac- 
terized by an exponent in between that of standard and 
single-file diffusion. Explicitly 



Sit) ~ A^ 



21+/J 



(i + ;3)r(V) 

which agrees with [241 ES] up to scaling. 



(i/T)(^+i)/' oc i(^+i)/2 (19) 



particle in a single-file diffusion system. We used a new 
simple method based on a quasi-equilibrium assumption 
and predicted the scaling behaviour and crossover times 
known from literature. A future challenge will be to try 
to improve our method to also yield correct prefactors. 
The reason to why the approach works so well is that 
particle density relaxations decays faster than the mo- 
tion of a tracer particle. It would be interesting to see 
how this assumption holds in other settings such a disor- 
dered single file where the diffusion constants are assigned 
randomly. 
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